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ABSTRACT 

Approximately 30 — 40% of all baryons in the present day universe reside in a warm-hot intergalactic 
medium (WHIM), with temperatures between 10^ < T < 10^ K. This is a generic prediction from 
six hydrodynamic simulations of currently favored structure formation models having a wide variety of 
numerical methods, input physics, volumes, and spatial resolutions. Most of these warm-hot baryons 
reside in diffuse large-scale structures with a median overdensity around 10 — 30, not in virialized objects 
such as galaxy groups or galactic halos. The evolution of the WHIM is primarily driven by shock 
heating from gravitational perturbations breaking on mildly nonlinear, non-equilibrium structures such 
as filaments. Supernova feedback energy and radiative cooling play lesser roles in its evolution. WHIM 
gas is consistent with observations of the 0.25 keV X-ray background without being significantly heated 
by non-gravitational processes because the emitting gas is very diffuse. Our results confirm and extend 
previous work by Gen & Ostrikcr and Dave et al. 

Subject headings: Gosmology: observations, large scale structure of Universe, intergalactic medium 
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1. INTRODUCTION 

Observations indicate that most of the baryonic matter 
in the universe does not reside in galaxies. At high red- 
shifts (z ^ 2), the overwhelming majority of baryons are 
in a diffuse, photoionized intergalactic medium (IGM), ob- 
servable as H I absorpt i on lines in the spectra of distant 



is the mean baryonic density) and temperature T: 



quasars (|Cen et al. 1994 ; Zhang. Anninos, fc Norman 1995 



Miralda-Escude et al. 1996| ; [Hernquist et al. 1996|) The 



baryonic d ensity inferred from these observations ( Rauclj 



et al. 1997) is in good agreement with nucleosynt hesis ar 



guments based on ob served deuterium abundances ( Tytler 



Fan, & Buries 1996). But by redshift zero the total bary- 
onic component, inferred from H I absorption, gas and 
stars in galaxies, and other observations, has declined to a 
level small compared to that see n at high redshift and ex- 
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Diffuse: 6 < 1000, T < 10^ K. Photoionized inter- 
galactic gas that gives rise to Lyman alpha absorp- 
tion. 



Condensed: 
galactic gas. 



S > 1000, T < 10^ K. Stars and cool 



Hot: T > 
groups. 



10^ K. Gas in galaxy clusters and large 



Warm-Hot: 10^ < T < 10^ K. The "Warm-Hot In- 
tergalactic Medium" (WHIM), discussed here. 



pectations from nucleos ynthesis ( Fukugita, Hogan & Pee 



bles 199q ; Hogan 1999). Thus the question arises, where 



Gen & Ostrike r (1999: hereafter |C099|) and Dave et al. 
(1999; hereafter DHKW ) predicted that a sizeable frac- 
tion of all baryons at the present epoch reside in this last 
warm-h ot phase (see Figure 2 of G09£ and Figure 12 of 



DHKW). Such a reservoir has significant implications for 



producing an accurate census of bar yons for comparison 
with nucleosy nthesis arguments (e.g. Fukugita, Hogan & 



Peebles 1998[) because gas at these temperatures and den- 



are the baryons at the present epoch? 

By the current epoch, hierarchical structure forma- 
tion has produced deep potential wells into which the 
baryons accrete, thereby moving a significant portion of 
the baryons from the IGM into stars, galaxies, groups, and 
clusters. These complex evolutionary processes may now 
be modeled directly using cosmological hydrodynamic sim- 
ulations, enabling an investigation into the location and 
phase of baryonic constituents in the present-day universe, 
and suggesting possible avenues for their direct detection. 

These hydrodynamic simulations of structure formation 
indicate that baryons in the universe reside in four broad 
phases, defined by their overdensity S = p/ p ~ 1 (where p 

^Princeton University Observatory, Princeton, NJ 08544; rad,cen,jpo@astro. princeton.edu 
^Spitzer Fellow 

^Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139; gbryan' 
"Hubble Fellow 

^Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138; lars@cfa.harvard.edu 
^Department of Astronomy, University of Massachusetts, Amherst, MA, 01003; nsk@kaka.phast.umass.edu 
'^Department of Astronomy, Ohio State University, Columbus, OH, 43210, dhw@astronomy.ohio-state.edu 
^Astronomy Department, University of Illinois at Urbana-Champaign, Urbana, IL 61801; norman,bwoshea( 

1 



sities is difficult to detect in either absorption or emission, 
as we will discuss in §^ 

In this paper we study the nature and evolution of 
warm-hot gas in a AGDM universe, using a suite of cos- 
mological hydrodynamic simulations having a wide range 
of numerical and physical parameters. The purpose of this 
paper is to ask how robust simulation predictions are to 
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Table 1 
Simulation parameters. 





Code 


n 




n 




Ho 




L 


eb 




Physics 


^WHIM 
Ob 


CwHIM ° 


Dl 


PTreeSPH 


0.4 


0.6 


0.95 


0.0473 


65 


0.8 


50 


7 


8.5 X 10® 


1,3,4 


0.30 


244 


D2 


PTreeSPH 


0.4 


0.6 


0.95 


0.0473 


65 


0.8 


11.11 


3.5 


1.1 X 10* 


1,3,4 


0.29 


405 


CI 


TVD-PM 


0.37 


0.63 


0.95 


0.049 


70 


0.8 


100 


200 


1.6 X 10® 


1,2,3,4 


0.42 


34 


C2 


TVD-PM 


0.37 


0.63 


0.95 


0.049 


70 


0.8 


50 


100 


2.0 X 10^ 


1,2,3,4 


0.37 


106 


Bl 


AMR 


0.3 


0.7 


1.0 


0.04 


67 


0.9 


100 


50 


7.9 X 10^ 


adiabatic 


0.32 


208 


B2 


AMR 


0.3 


0.7 


1.0 


0.04 


67 


0.9 


100 


1 


9.9 X 10® 


1,2,3,4 


O.St 


400t 



''Box size in eomoving h ^Mpc. 

''Spatial resolution in eomoving /i^'kpe; for PTreeSPH and AMR, this is the highest resolution achieved in dense regions. 

'^Baryonic mass resolution in Mq. In TVD-PM, this is the average mass per cell. 

''l = H, He cooling; 2 = Metal cooling; 3 = Photoionization; 4 = Star formation & feedback. 

"Clumping factor of warm-hot gas at 2 = 0; see ^ 

^Values extrapolated from z = 0.75 to z = based on a comparison with simulation Bl. 



these parameters, to examine the physical state of warm- 



hot gas in the universe, and to investigate constraints on 
warm- hot gas from soft X-ray background observations. 

In §H we briefly describe the simulations used. The evo- 
lution of diffuse gas at high redshift into condensed, hot, 
and WHIM gas at the present epoch is quantified in §|[ 
Our primary result, presented in §^ is that the WHIM ac- 
counts for a significant fraction 30 — 40%) of baryonic 
mass at z = 0, regardless of variations in spatial resolution, 
input physics, or hydrodynamic algorithm. This is because 
the evolution of the WHIM is driven primarily by shock- 
heating of gas falling into gravitationally-generated poten- 
tial wells, a process which is well-understood and mod- 
eled, and only secondarily by processes such as supernova 
feedback, radiative cooling, and photoionization. Further- 
more, most of the WHIM gas is at relatively low overden- 
sities, so shock heating of intergalactic gas occurs during 
flows onto non-equilibrium large-scale structures such as 
filaments. The majority of warm-hot gas is found out- 
side of virialized objects such as galactic halos and galaxy 
groups. The low overdensities explain why radiative cool- 
ing and supernova heating do not drive its evolution, and 
why the presence of this component probably does not vi- 
olate constraints from the X-ray background, as we show 
in Finally, we briefly discuss strategies for direct de- 
tection of this gas, noting that the easiest place to de- 
tect emission from WHIM gas is relatively close to galax- 
ies, where it is dense, even though most of the warm-hot 
baryons are not in these regions. In summary, the WHIM 
is a robust and generic prediction of currently popular 
hierarchical structure formation models, and it contains 
roughly one-third of all baryons in the universe today. 

2. SIMULATIONS 

We use six cosmological hydrodynamic simulations of 
randomly-selected volumes in A-dominated Cold Dark 
Matter universes, employing three different numerical 
techniques, with a range of physical and numerical pa- 
rameters. These parameters are summ arized in Table 1. 
Simulations run with Parallel TreeSPH (Dave, Dubinski fc 



Hernquist 1997 ) are labeled Dl and D2, simulations run 
with TVD-PM ( |Ryu et al. 199^ ) are labeled CI and C2, 



and simulat ions using Adaptive Mesh Refinement (AMR; 
Bryan 1999| ) are labeled Bl and B2. The cosmology chosen 
is close to the "concordance model" which is in agreeme nt 
with a wide variety of observations ( Bahcall et al. 1999| ). 

For our purposes, the significant inputs are the vari- 
ations in spatial resolution (1 200/i~^kpc), baryonic 
mass resolution (2 x 10^ lO^M©), boxsize (11.11 
100 /i^^Mpc), input physics, and hydrodynamic algo- 
rithms. Dl and D2 are high spatial resolution Lagrangian 
(particle-based) simulations, CI and C2 are lower spa- 
tial resolution Eulcrian simulations having high mass res- 
olution and employing the Total Variation Diminishing 
scheme, and Bl and B2 are high resolution adaptive mesh 
simulations based on the Piecewise Parabolic Method 
(note the exceptional resolution of ~ 1/i^^kpc achieved by 
AMR in the high density regions). All simulations except 
Bl include radiative cooling from H and He, photoion- 
ization heating, and star formation; Bl includes none of 
these. CI, C2, and B2 additionally include metal-line cool- 
ing, with the metallicity determined self-consistently from 
supernova output. All simulations use a A-dominated cold 
dark matter universe, having similar power at cluster and 
galaxy scales. B2 has only been evolved to z = 0.75, but 
this will be sufficient to indicate the relevant trends. 



3. EVOLUTION OF INTERGALACTIC GAS 

Figure 1 shows the evolution of the baryonic mass frac- 
tion in the four phases described above. The four panels 
show results from simulations Dl, D2, CI and C2. Despite 
differences in simulation volume, resolution, and numeri- 
cal method, the evolution of various phases is qualitatively 
similar. At high redshift (z > 2), the dominant fraction 
of baryons resides in diffuse gas (dashed lines) giving rise 
to Lyman alpha fore st absorbers, as has been explored in 
detail elsewhere (see Ranch 199^ for a review). As struc- 
ture forms, diffuse gas is shock-heated, producing warm- 
hot gas (solid lines). Gas that is driven to higher densi- 
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ties due to gravitational instability is able to cool into the 
condensed phase (dotted hnes) and form stars. At lower 
redshifts {z <^ 2), large potential wells are produced that 
shock-heat gas to T > 10'' K, giving rise to hot cluster gas 
(dot-dashed curves). Nevertheless, at the present epoch, 
the total fraction of baryons in clusters is small. 
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Figure 1: Evolution of mass fractions in four baryonic phases, 
in four simulations. 

Figure 1 shows that the fraction of diffuse baryons at 
the present epoch is between 20% and 40%, the fraction 
of warm-hot baryons is 30 — 40%. Gravitationally bound 
baryons, i.e. those in stars, galactic gas, clusters, and in- 
tragroup media, make up the rest. Thus baryons in the 
present-day universe are divided roughly equally among 
diffuse, warm-hot, and bound components. 

While all the simulations are qualitatively similar, the 
exact distribution of gas in these phases is sensitive to de- 
tails of the simulations. For instance, the condensed phase 
fraction is sensitive to how gas cools and forms stars in 
these simulations, which in turn is significantly affected by 
resolution (since cooling and star formation oc p^). Also, 
the effect of supernova feedback heating depends strongly 
on resolution, as we will discuss later. The growth of struc- 
ture is affected by the amount of large-scale power present. 
In particular, the hot gas fraction in clusters is sensitive 
to cosmic variance, since it is dominated by the largest 
virialized objects in the volume. Other phases, including 
the warm-hot phas e, are less sensitive to volume effects, 
as we show in 



4.1 



All things considered, it is not surprising that there are 
differences up to a factor of two in the fractions in various 
phases at z = 0. While the differences may be signifi- 
cant, an investigation of their exact causes is beyond the 
scope of this work, though these differences offer clues into 
the physica l pr ocesses driving WHIM evolution, as we will 
explore in §0. Rather, we focus on the remarkable qual- 
itative consistency in the evolution of various gas phases, 
given the variety of simulation methodologies utilized. 

For the condensed phase, a census of baryons in stars 



and cold gas estimates its mass fract ion to be around 
20% dFukugita, Hogan fc Peebles 19981 ). The simulations 
shown generally produce somewhat higher values for the 
condensed fraction, but there are uncertainties in the con- 
tributions from low-mass stars and supernova remnants. 
Of these simulations, C2 contains the largest fraction of 
condensed gas (~ 40%), while the rest are all lower, down 
to ~ 25% for simulation Cl. Bl, of course, has virtually 
no condensed gas since it does not include cooling. In 
what follows, we will fix the condensed fraction to be 20% 
and redistribute the excess condensed gas equally among 
all the other components, assigning a somewhat arbitrary 
fraction of 1/3 of it to the WHIM (the only component 
we examine from here on), in order to facilitate a less 
resolution-dependent comparison of the intergalactic com- 
ponents in these simulations. 

4. THE WARM-HOT INTERGALACTIC MEDIUM 
4.1. Evolution of flwHiM 




Figure 2: Evolution of Hwhim in our six simulations. For this 
comparison, the baryonic fraction in cold galactic gas and stars 
has been fixed at 20% in all simulations. 

Figure 2 shows the evolution of the mass fraction of 
baryons in warm-hot gas, ilwHiM/^^b, for our six simula- 
tions, with the condensed component having been fixed 
at 20%. Cl and C2 have the highest WHIM fractions, at 
42% and 35% of baryonic mass, respectively, at z = 0. Dl 
and D2 have lower WHIM fractions at the present epoch, 
around ^ 30%. B2 has only been evolved to z = 0.75, but 
its evolution closely mirrors that of the Lagrangian runs, 
so we expect its WHIM fraction will also be ^ 30% by 
z — 0. These values are listed in Table 1. By the present 
epoch, all simulations have WHIM fractions within 50% of 
each other. Simulation Bl does not include radiative cool- 
ing, so it is an unphysical model that will only be used for 
comparison with B2. 

We expect that the evolution of warm-hot gas is gov- 
erned by shock heating of intergalactic gas onto large-scale 
structure, supernova feedback, and radiative cooling. Nu- 
merical considerations such as resolution, volume and hy- 
drodynamic algorithm may also play a role. In this sec- 
tion we use case-by-case comparisons among our six sim- 
ulations to examine how each of these processes affects 
WHIM evolution in our models. 

All of our simulations explicitly include the growth of 
structures and the accretion of gas onto those structures. 
The qualitative consistency of WHIM evolution in all our 
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simulations, including simulation Bl without cooling or 
star formation, suggests that gravitational shock heating 
of gas falling on large-scale structures is the dominant 
heating mechanism for WHIM gas. This is an important 
point, as it suggests that all other differences between these 
simulations are of secondary importance. 

The most obvious differences between models are that 
Dl, D2 and B2 all predict similar WHIM evolution, while 
CI and C2 have more WHIM gas. The primary distinc- 
tion between these sets of simulations is spatial resolu- 
tion. This manifests itself in various ways, as it affects 
gas cooling, the rate at which stars form, and the the 
injection of supernova heat energy into intergalactic gas. 
To elaborate on this last point, high-resolution runs (Dl, 
D2, B2) deposit feedback energy (thermal only) locally 
in very high density regions where stars are forming, and 
thus it quickly radiates away, with virtually none of it be- 
ing distributed into the intergalactic medium. Conversely, 
lower-resolution runs (CI and C2) deposit the same su- 
pernova energy over hundreds of kiloparsecs, resulting in 
a significant fraction of feedback energy escaping into dif- 
fuse regions where it cannot radiatively cool away. None 
of the simulations here are capable of resolving supernova- 
driven galactic winds through a multi-phase interstellar 
medium, which are likely to be responsible for distributing 
energy and meta ls into the diffuse IGM (Mac Low 2000; 
Efstathiou 2000), thus we are relying on heuristic model 



ing of these processes. As we discuss in §4.2, it is in the 



distribution of feedback energy where resolution plays its 
most crucial role in predicting the evolution of WHIM gas. 

First we examine the effect of radiative cooling. Simu- 
lation Bl has no cooling, whereas simulation B2 is a simi- 
lar run with cooling (and is our highest spatial resolution 
run). Cooling has a greater effect at earlier times because 
intergalactic gas is denser then and thus can cool signif- 
icantly. After 2; 3, the rate of growth of the WHIM 
fraction is similar in Bl and B2, suggesting that WHIM 
gas is no longer affected by cooling. Even by z = 0.75, 
the difference between Bl and B2 is not large, indicating 
that cooling plays a minor role in the overall evolution of 
WHIM gas. 

Eulerian codes resolve shock fronts better than SPH in 



diffuse regions (Kang et al. 1994), because they have a 
higher density of resolution elements there, and, in the case 
of the Eulerian codes discussed here, because they incorpo- 
rate explicit shock capturing algorithms to resolve fronts 
over two cells. If a significant component of shock heating 
arises from small-scale shocks unresolved by SPH, Eule- 
rian codes may produce higher temperatures. One might 
suspect that the difference between PTreeSPH runs and 
TVD-PM runs could be partially due to this effect. How- 
ever, B2 is also an Eulerian simulation, yet it has nearly 
the same WHIM fraction as the Lagrangian runs. So this 
cannot be a significant effect. 

Simulation volume could also play a role in the amount 
of WHIM gas, since larger volumes contain larger per- 
turbations that can result in stronger gravitational shock 
heating. This may explain the difference between CI and 
C2, whose box lengths differ by a factor of two. However, 
Dl has a box length five times that of D2, yet their WHIM 
fractions are nearly identical at all times. Furthermore, B2 
has twice the box length of Dl, yet its WHIM fraction is 
in good agreement with Dl and D2. As we will show in 



the next section, typical WHIM gas is at very moderate 
overdensities, and cosmic variance in that regime is typi- 
cally small, in contrast to the hot IGM fraction which is 
dominated by rare, massive objects (i.e. clusters) whose 
numbers are more sensitive to simulation volume. 

It is possible that opposing effects are causing Dl, D2 
and B2 to be similar. For instance, Dl resolves shock 
fronts better than D2, perhaps resulting in more shock 
heating and making up for its lack of large-scale power. 
It is also possible that the adaptive refinement of B2 re- 
sults in more cooling in shock fronts, compensating for its 
larger volume. However, unless all of these effects are of 
secondary importance, cancellations at the level we find 
here would require a remarkable coincidence. 

The parameters of the underlying cosmological model 
are expected to have a non-negligible effect on the WHIM 
component. The rate of structure evolution in large part 
determines how much shock-heated gas is present at any 
epoc h. While we have not sampled a range of cosmologies 
here, DHKW examined four cosmologies, namely ACDM, 
Tilted CDM, Cold-|-Hot DM, and Open CDM models, us- 
ing numerical parameters similar to those of simulation 
D2. The differen ces betw een the WHIM fractions (called 
"shocked" gas in DHKW) in those models at the current 
epoch is comparable to the differences seen here due to 
other factors, as indicated by their Figure 12. Similarly, 
C099| examined ACDM, Cold+Hot DM, and Open CDM 



models, and also found broad consistency. Thus while our 
quoted fractions may be specifically for a ACDM cosmol- 
ogy, our qualitative conclusions are unlikely to be highly 
sensitive to cosmology. 

In summary, the evolution of rJwHiM is qualitatively 
consistent among all simulations examined, and results in 
w 30 — 40% of baryons residing in the WHIM today. Ra- 
diative cooling, simulation volume, and algorithmic details 
do not significantly affect WHIM evolution in these runs. 




Figure 3: WHIM gas in simulation C2. Contours are color 
coded by overdensity; green represents overdensity S ~ 10, 
while red shows S ~ 10*. 
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4.2. The Physics of WHIM Gas 

We now explore the physics that drives the formation 
and evolution of WHIM gas. A qualitative physical pic- 
ture of WHIM gas may be obtained by examining Fig- 
ure 3, which shows the location of gas in simulation C2 
having 10^ < T < lO^K, with contours color-coded by 
density. WHIM gas is seen to primarily trace out filamen- 
tary large-scale structures. Like the intergalactic medium, 
WHIM gas does cluster around dense regions that are sites 
of galaxy formation. However, we will show below that the 
majority of WHIM gas is contained in the filaments. 
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Figure 4: Top panel shows the mass fraction of WHIM gas as a 
function of density at z = in our simulations. Bottom panel 
shows the same quantity a,t z — 0, 1, 2, 3 in simulation Dl. 

The top panel of Figure 4 quantifies the spatial distribu- 
tion of WHIM gas in the universe. It shows a histogram of 
WHIM gas mass as a function of density, for simulations 
Dl, D2, CI, and C2. All four simulations consistently 
show that the dominant fraction of WHIM gas is at rela- 
tively low densities, with a peak around an overdensity of 
~ 10 — 30. Simulations Bl and B2 are not shown because 
Bl does not include cooling and B2 has only been evolved 
to z = 0.75. StiU, Bl at z = and B2 at z = 0.75 show 
peak overdensities of 18 and 21, respectively, so they are 
consistent with the other simulations. 70 — 80% of WHIM 
baryons lie in the overdensity range 5 < 5 < 200, typical 
of filaments. 

The typical overdensities of WHIM gas are much smaller 
than that of matter contained in bound, virialized objects. 
A maximal estimate of the bound fraction of WHIM gas 
may be obtained as the fraction of WHIM gas with (5 ^ 60, 
which is approximately the overdensity at the virial radius 
of an isothermal sphere in a ACDM cosmology. The mass 
fraction of WHIM gas with 5 ^ 60 is « 30% in our sim- 
ulations. Clearly some of this gas will not be bound, but 



rather infalling material. As an independent check, we use 
the group finding algorithm SKID^ (Spline Kernel Inter- 
polative DENMAX) to identify bound warm-hot particles 
in simulation Dl, and find a bound fraction of WHIM 
gas between ~ 10% and 25%, depending on the linking 
length used (50 — 500/i~^kpc). The lower end of this range 
probably corresponds to gas contained in galactic halos, 
although the extent of a galactic halo becomes ill-defined 
when it resides within a group or cluster. In summary, 
WHIM gas is mostly an intergalactic component, with a 
majority of it residing outside of virialized structures such 
as galaxies or groups. Note that coronal gas in galaxies 
also lies in the warm-hot temperature range, but it is a 
very small fraction of the galactic baryonic mass, as most 
galactic baryons are tied up in stars and cold gas. 

How can intergalactic gas be heated to T > 10^ K with- 
out being bound in a massive virialized halo? One way 
would be if energy was added from non-thermal processes 
such as supernova feedback. However, this is not the driv- 
ing process for WHIM gas in these simulations. Supernova 
feedback energy is added in all simulations shown in Fig- 
ure 4, but in simulations Dl, D2 and B2 the effect is negli- 
gible, since, as discussed before, feedback is added purely 
thermally into very dense regions, where it radiates away 
almost immediately and adds no heat to diffuse regions. 
But even in these high-resolution simulations, where su- 
pernovae add negligible heat to WHIM gas, the typical 
overdensities are smaller than those typical of virialized 
halos. 

Supernova heating is, however, likely to be responsible 
for the somewhat lower overdensities in simulations CI 
and C2 as compared to the higher resolution simulations, 
so its effect is non- negligible. The differences cannot be 
due to other spatial resolution effects such as cooling, as 
CI and C2 are themselves quite similar, as are Dl and D2. 
Instead, the lower spatial resolution in CI and C2 results 
in significant supernova energy being deposited in inter- 
galactic g described in the previous section. This 
additional feedback heating raises the pressure of inter- 
galactic gas and lowers the typical density. Furthermore, 
because supernova energy is not radiated away immedi- 
ately, there is a small component of warm-hot gas in and 
around galaxies seen as the high-density bump for CI and 
C2 in Figure 4. Note that the densities in Dl and D2 are 
computed slightly differently than in CI and C2; in the 
former, Lagrangian runs, the density field is smoothed by 
the (variable) smoothing length of the particle, while in the 
Eulerian runs the smoothing is fixed at the cell size. How- 
ever, since in both cases the smoothing is done on scales 
smaller than the density variations, this is not expected to 
produce any systematic differences in the resulting densi- 
ties. 

The effect of supernova heating can be roughly esti- 
mated by the following simple argument. Consider two 
extreme cases, one where all supernova energy is deposited 
within dense regions where it immediately radiates away, 
and another where supernova energy is distributed uni- 
formly over all baryons in the universe. The former case 
is a reasonable approximation for the Lagrangian runs Dl 
and D2 as well as the high-resolution AMR run B2, while 
the latter is closer to the Eulerian runs CI and C2 (though 
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clearly a much more extreme case). In the former case, su- 
pernovae add no heat to WHIM gas. In the latter, the spe- 
cific heat added per baryon yields a temperature increase 
of 



SN 



^* Pcrit eSN 



(1) 



0.03 



0.02 - 



0.01 




where fi* is the cosmic mass fraction in stars, pcrit is the 
critical density, esN is the specific energy output per unit 
mass of stars formed, and hh is the mean number den- 
sity of H atoms. We take f2* = O.lflb = 0.002ft,~^, and 
esN = 2.5 X 10^^ erg g'^Mg from a Salpeter IMF with 
each supernova from a star with M > 8Mq outputting 
10^1 ergs. Then, STsn = 0.2 keV « 2 x 10^ K. Assuming 
isentropic heat distribution, pT^ = constant, so gas with 
T « 4 X 10^ K (the maximum of the temperature distri- 
bution, as we will show in Figure 5) will have its density 
reduced by roughly a factor of two compared to the no 
heating case. This is roughly the level of reduction seen in 
the Eulerian as compared to the Lagrangian runs. Clearly 
this model is overly simplistic, as it would predict that 
all intergalactic gas has T ^ 10^ K, but it does roughly 
indicate the magnitude of the effect of supernova heating. 

Note that there are no other heating processes in these 
simulations that contribute significantly to the WHIM. For 
instance, photoionization only heats gas to ~ 10^ K. Any 
other heating process one could devise, such as cosmic 
rays or supernovae occuring in dwarf galaxies in voids, 
would add more pressure support to the gas, and thus 
would push the typical density of WHIM gas even lower. 
Thus our simulations indicate that WHIM gas is heated 
to T > 10^ K primarily by shock-heating of gas accret- 
ing onto large-scale structure. These structures, typically 
filamentary, are not virialized or in dynamical equilibrium. 

In the previous section we showed empirically that ra- 
diative cooling has a minor effect on WHIM evolution. 
This may be understood physically given that warm-hot 
gas is typically at such low overdensities that it cannot col- 
lapse into virialized objects, and due to the metagalactic 
photoionizing background most inter galactic gas lies well 
within the optically thin regime (see DHKW, Figure 10). 
Thus typical WHIM gas is too diffuse to self-gravitate or 
self-shield, and has no way to achieve the densities required 
to make radiative processes important. Improving the nu- 
merical resolution of our simulations would not change this 
result, as it is based on simple physical arguments. 

The bottom panel of Figure 4 shows a histogram of 
WHIM gas mass for the Dl model at z = 0, 1, 2, 3 (solid, 
dotted, dashed, long dashed hues). There is a shght trend 
to higher overdensities at earlier times, since at early times 
the contribution from virialized structures is greater. This 
is because gas has not had time to accrete and shock on 
more diffuse structures, and the largest structures at ear- 
lier times have lower temperatures that can fall into the 
warm-hot range. Still, this is a minor effect; basically, the 
peak overdensity does not evolve significantly with red- 
shift. 



Figure 5: The mass fraction of baryons as a function of tem- 
perature in simulation Dl, at z = (solid), z — 1 (dotted), 
z = 2 (short dashed) and z — 3 (long dashed). The arrows in- 
dicate the predicted peak temperature from gravitational shock 
heating at various z, from equation ^. 

Figure 5 shows the temperature distribution of inter- 
galactic gas with T > 10^-^ K in simulation Dl at z = 
0,1,2,3. The amount of gas at these warm- hot and hot 
temperatures grows with time, as seen from Figures 1 and 
2. The temperature at the peak of the distribution grows 
in time as well, reflecting the fact that the universe con- 
tains hotter structures at later times. The peak tempera- 
ture at all redshifts falls within the warm-hot range, and 
by z = it is up to ~ 4 X 10^ K. Note that this is close 
to the temperature of th e excess diffuse emission seen by 
Wang fc McCray (1993)| in ROSAT data. Figure 5 also 
shows that the mass of the WHIM component is insen- 
sitive to our somewhat arbitrary choice of 10^ K as the 
defining lower temperature. A choice of 10^'^ K (as in 
DHKW) or even 10^'^ K would not drastically affect our 



conclusions. 

An analytic estimate of the evolution of the peak tem- 
perature can be obtained by considering the temperature 
of intergalactic gas shock heated on mildly nonlinear large- 
scale structure. If the length scale going nonlinear at a 
given epoch is L^i, and the perturbation collapses on a 
timescale t, the resulting sound speed behind the shock 
will be ~ Lni/t. Given that the perturbation has taken a 
Hubble time to collapse, t ~ H~^, where H is the Hubble 
constant at time t. The resulting post-shock temperature 



is then (0099) 



(2) 



where K is a constant, at any given epoch. The value of 
Tni, with K = 0.3, is shown at each redshift plotted in 
Figure 5 by the arrows above each curve. This value of K 
produces roughly the correct peak temperature at z = 0, 
and as can be seen from Figure 5, it is a reasonable fit 
to the peak temperature at higher redshifts, though it be- 
comes progressively more difficult to identify a peak. This 
value of K also produces the correct evolution of the glob- 
ally aver aged te mperature in simulations CI and C2, as 
shown in C099. Thus the evolution of the temperature 



distribution of WHIM gas is consistent with the interpre- 
tation that it is heated by gravitationally-induced shocks 
on mildly nonlinear large-scale structure. Since there is a 
wide range in the properties of the collapsing structures 
and therefore infall velocities, there is also a wide range in 
gas temperatures. 
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Figure 6: Contours in temperature and density for simulation 
Dl at z = 0, enclosing 10%, 50% and 90% of the baryons in 
the range shown. Density and temperature are correlated in the 
WHIM regime. Thick line indicates a scaling oi p/ pb — T /l^""^ 
in the warm-hot temperature range. 

The temperature and density of WHIM gas are corre- 
lated. Figure 6 shows a contour plot of mass within the 
warm-hot range, as a function of density and tempera- 
ture, for the Dl model at z = 0. The contour levels en- 
close 10%, 50% and 90% of the mass in the temperature 
and density ranges shown in the plot. The thick hne in- 
dicates an "equation of state" p oc T that provides a rea- 
sonable fit to gas in the range 10^ < T < lO'^ K. This 
relationship is different from that of diffuse gas, which 
typically has p oc T^'^, and the temperature-density re- 
lation of WHIM gas has much greater scatter. The higher 
temperature, different slope, and greater scatter all reflect 
the importance of shock heating as the dominant mech- 
anism controlling the thermal properties of WHIM gas; 
the "equation of state" for diffuse gas, on the other hand, 
arises from the competition between photo ionization and 
adiabatic cooling due to Hubble expansion ( Hui & Gnedin 
19970 Figure 6 also suggests that detecting WHIM gas in 
emission will be easier for gas that is at the highest end 
of the WHIM temperature range, since that gas will be 
both denser and hotter. However, the dominant portion 
of WHIM gas is at lower temperatures, which is perhaps 
most easily de tected via absorption lines (Tripp, Savage & 
Jenki js 200"o|) . 



5. CONSTRAINTS FROM THE SOFT X-RAY BACKGROUND 

Gas with temperatures in the range 10^ < T < 10^ K 
will emit thermally in the soft X-ray band. The extra- 
galactic soft X-ray foreground (SXRB) flux at 0. 1-0.4 keV 



is roughly ^ 20 — 35 keV cm s sr keV ( Warwick 



& Roberts 1998), though uncertainties are large because 
galactic coronal gas provides an increasing foreground to 
lower energies. At slightly higher energies (~ 1 keV), 
the XRB has further been resolved nearly co mpletely 
80 - 90%) i nto point sources, mostly AGN ( [Mushot-] 
sky et al. 2000). Reasonable arguments then allow only 



a small contribution to the SXRB from diffuse gas, ^ 
4 keV cm-V_^r~ikeV"i (Wu, Fabian & Nusser 1999, 



hereafter WFN). Such a limit, in principle, places con- 
strain ts on the amount of gas in the universe at warm-hot 



tempc ratures. 



These limits were explore d in tw o independent papers 
using similar methodologies, WFN and Pen (1999). Both 
papers argue that the standard picture of hierarchical for- 
mation of viriahzed objects results in a predicted SXRB 



that exceeds the observed limits. They suggest that sig- 
nificant non-gravitational heating, typically ~ 1 keV per 
baryon, is required to unbind warm-hot gas from virial- 
ized objects in order to satisfy the SXRB constraints. In 
this section we discuss these constraints in the context of 
WHIM gas, and find that our simulations paint a very 
different picture for soft X-ray emission than the simple 
models assumed in those two papers. A full calculation 
of the SXRB from these simulations is a complicated un- 
dertaking (because of metallicity, bandpass, and numerical 
issues) that is beyond the scope of this paper. However, 
the physical properties of WHIM gas in the simulations 
are quite different from the properties assumed by WFN 
and Pen (1998), and we will show that scaling our results 
to theirs suggests that the observed SXRB fiux does not 
rule out the WHIM predicted by these simulations. The 
essential difference is one of density . Both WFN a nd Pen 
(1998) base their calculations on a Press-Schecter (1974) 
analysis, which implicitly assumes that gas is in virialized 
objects with typical overdensity ^ 200. However, most of 
the WHIM gas in the simulations is in lower-density fila- 
mentary structures rather than virialized objects, thus the 
SXRB emission is lower (cf. Figure 5). 

We can quantify this difference in typical overdensity by 
considering the clumping factor of the emitting gas. If we 
define the clumping factor for gas component g as 



Co =< Pg > / < Pg >^ 



^9 =^ l^a ^ / ^ h'g ^ J (3) 
then the free-free emissivity from that component is 

esxRB (X < p^tO-S > cx Cgfi^TO ■^ (4) 

where pg, Slg, Tg « 10^ K, and Cg are the density, mass 
fraction, temperature, and clumping factor of the gas emit- 
ting in soft X-rays. The fiux j'sxrb of soft X-ray back- 
ground is t hen th e emissitivity multiplied by path length 



WFN argue, sensibly, that it is predominantly 
warm-hot gas (10^ < T < 10^ K) that is responsible for 
soft X-ray emission. This means the appropriate clumping 
factor Cg is that of warm-hot gas, Cwhim- 

There are several ways to calculate Cwhim in our sim- 
ulations. One can directly calculate it from equation ^ 
which is the approach we use for our Eulerian simulations. 
For Lagrangian simulations, because each particle repre- 
sents a different volume of gas, it becomes more numeri- 
cally convenient to calculate Cwhim ~ Cwhim (0), where 
Cwhim (?") is the two-point correlation function of WHIM 
gas at radius r. 

The resulting Cwhim values for our simulations at z = 
are listed in Table 1. All simulations show clumping fac- 
tors in the range ~ 30 — 400. The smaller clumping factors 
in CI and C2 arise directly because the WHIM gas is typi- 
cally less dense in these models as compared to Dl, D2, Bl 
and B2 due to greater supernova feedback in the diffuse 
IGM, as explained before. 

Our cl umping fact ors are significantly less than that 



used by [Pen (1999)| , who adopts C ^ 900 and argues 
10**. The reason is that the clumping factor of Pen 



C > 



(1999) is actually that of all gas (assumed to trace the dark 



matter), not the warm-hot baryons. Such a clumping fac- 
tor is dominated by the contribution from galaxies and 
other collapsed, virialized structures. This is not the ap- 
propriate clumping factor with which to calculate the soft 
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X-ray background, since as we have shown in the previous 
section, WHIM gas is typicaUy at much low er overdensi- 
ties than viriahzed structures. In Pen (1999) , the clumping 
factor is strongly dependent on simulation resolution, be- 
cause it is dominated by the highest overdensity objects. 
Amongst our simulations, a wide range of resolutions give 
barely an order of magnitude difference in clumping factor, 
since the WH IM gas is predominantly diffuse. 

FcnXTUUU) determines that ^ 1 keV of non-gravitational 
heating would reduce his clumping factor to C ^ 60, which 
he argues is necessary to satisfy XRB constraints. How- 
ever, this was derived assuming that all gas contributes to 
the SXRB (see his eq. 2), whereas in our simulations only 
~ 30% of the baryons contribute to the SXRB (probably 
an even smaller fraction since the narrow range of SXRB 
energies arises from a narrower range of gas temperatures 
than 10^ < T < lO'^ K). From equation |[ we see that a 
reduction in fig, the mean density of emitting gas, results 
in an increase in the maximum allowed clumping factor 
by (ilfc/OwHiM)^ ~ 10. Thus his analytically-derived con- 
straint translates to Cwhim <^ 600. Our clumping factors 
(< 400) are consistent with his analysis, despite the fact 
that WHIM gas in our simulations undergoes almost no 
non-grav itational heating. 

WTW present a more detailed model from which we can 
crudely estimate our SXRB by comparison. They do not 
specifically use clumping factors to predict the SXRB in 
th eir model, b ut they perf orm a similar ca lculation to that 

analysis with an 
formalism makes 



The two-point correlation function ^whim illustrates 
quite clearly that WHIM gas is not bound in virialized 
halos, and is instead associated with large-scale structure. 
Figure 7 shows ■^whim for simulation Dl (solid line), D2 
(dashed), as well as f galaxy (dotted line) calculated from 
Dl- '^galaxy is computed as the correlation function of all 
star particles in this simulation, where stars are formed 
slowly out of cold, de nse gas and are seen to trace th e 
galaxy population (see Katz, Weinberg & Hernquist 1996). 



of |Pen (1999)1 , based on a jPress-Schcctcr 



assumed halo profile. The Press-Schcctci 
two assumptions about how gas is distributed and heated 

in the universe: fl) All soft X-ray emitting gas is bound in versus p revious studies arise from their use of a [Press 



At scales less than a few hundred kpc, ^whim flattens, in- 
dicating that WHIM gas does not cluster on scales smaller 
than those typical of large-scale structures. Conversely, 
the galaxy correlation function continues as approximately 
a power law in to the resolution limit of our simulation. 
The clumping factor of WHIM gas, ^whim(O), is well- 
defined since we can meaningfully extrapolate from our 
resolution limit to r = 0. This is not true for the clumping 
factor of galaxies, for which we may only set a lower limit 
(consistent with Pen 1999), Cgaiaxy > 4 x lO**. 

We can crudely estimate the SXRI3 flu x we w ould pre- 
dict in comparison to that predicted by [WFN by multi- 
plying their predicted jsxrb by CwHiM/Cgaiaxy This cor- 
rects for the fact that they use virialized objects (which we 
take as galaxies) instead of the true diffuse WHIM respon- 
sible for SXRB emission. Thus we would roughly predict a 
SXRB flux th at is tw o orders of magnitude lower than that 
predicted by WFN. Examining their Figure 2 or 3 shows 
that such a reduction in j'sxrb makes the SXRB contri- 
bution from warm-hot gas consistent with observational 
limits. 

In summary, the discrepancies in the predicted sof t XRB 



viriali bcd halos (namely, groups of galaxies): (2) Gas can Schecter] analysis, or an equivalent method that implic- 



only be heated by accretion onto a virialized halo, with the 
gas te mperature se t by the virial temperature of the halo. 
Unlike Pen (1999), they allow some gas to cool and there- 
fore not emit in X-rays, but this is a small correction. In 
their model, the gas clumping factor is the clumping factor 
of virialized objects, since they are assuming that emitting 
gas solely resides in such objects. 



Our simulations suggest that both [Press-Schectei| as- 
sumptions are strongly violated in the case of WHIM gas. 
First, only a small fraction of WHIM gas is bound in viri- 
alized halos; most is distributed much more diffusely. Sec- 
ond, in our simulations, intergalactic gas can be heated 
significantly by purely gravitational processes prior to be- 
ing accreted onto a virialized object. 




0.05 0.1 0.5 
r (h^'Mpc) 



Figure 7: Two-point correlation functions for WHIM gas from 
simulation Dl (solid line) and D2 (dashed line), and for galaxies 
from simulation Dl (dotted line). 



itly places emitting gas into virialized objects, in order to 
study an intergalactic component of baryons that is in- 
herently much more diffuse. These methods do not allow 
for gravitational shock heating in unbound objects such 
as filaments, thus they are forced to postulate significant 
non-gravitational heating to make the emitting gas more 
diffuse. In our scenario, most WHIM gas has never fallen 
into virialized objects, though it is still heated (almost) 
purely via gravitational processes. Scaling previous anal- 
yses by the typical overdensities found in our simulations 
suggests that the properties of the predicted WHIM are 
consistent with SXRB constraints. However, these scal- 
ing arguments are approximate at best, and not enough 
to guarantee consistency. More accurate calculations of 
X-ray emission from these simulations are certainly war- 
ranted and may yield interesting constraints on the WHIM 
component; these issues will be addressed in detail in sep- 
arate papers (Fardal et al., in preparation; Phillips et al., 
in preparation). 

6. SUMMARY 

We study the warm- hot intergalactic medium (WHIM), 
defined as all the gas in the universe with temperature 
10^ < T < 10^ K, in six cosmological hydrodynamic sim- 
ulations with widely varying spatial resolutions, volumes, 
code algorithms, and input physics. In each simulation, 
the WHIM contains « 30 - 40% of all baryons in the 
present-day universe. As a rule of thumb, our simula- 
tions predict that the fractions of baryons in the warm- 
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hot phase, the diffuse phase, and gravitationaUy bound 



systeris are roughly comparable at the present epoch. 

The WHIM is comprised primarily of gas at moder- 
ate overdensities, with a median overdensity fO — 30, 
and is predominantly an intergalactic component. It is 
shock-heated by accretion onto non-equilibrium filamen- 
tary large-scale structures, with possibly a small energy 
contribution from non-gravitational processes such as su- 
pernova feedback. Despite being heated primarily by grav- 
itational processes, WHIM gas in our simulations is consis- 
tent with constraints from the soft X-ray background. The 
clumping factors of WHIM gas in our simulations range 
from ~ 30 — 400, which are far below that of virialized 
objects 10*) that some studies assume are the sources 
of soft XRB photons. 

Different simulations give somewhat different fractions 



of O VI absorbers a t redshifts 0.14 <; z <; 0.27 ( [Tripp, Sav- 
age fc Jenkins 200C ; Tripp fc Savage 200C ) , which together 



with conservative ionization corrections and mctallicities 

implies riwHiM ^ ■^^'^-o'm2^75 y O'' ^wmu/^b ^ 10% at 
the present epoch. However, some of these O VI absorbers 
may be photoionized, arising in cooler, low-density inter- 
galactic gas at T 10* K. 

WHIM gas emission may be easiest to detect around 
high density regions such as clusters because that is where 
the density and temperature are highest within the WHIM 
range (cf. Figure 6), but this is not where the majority of 
WHIM gas is located. For instance, high-resolution Chan- 
dra spectra of regions between clusters that are free of 
bright sources may provide a detectable signal (Phillips 
et al., in preparation), and observations from XM M could 
possibly image filaments of warm-hot gas directly ( [Pierre,| 



of WI pIM gas at the present epoch, though the fraction is Bryan fc Gastaud 200Q|) . Our simulations, however, pre- 



significant in ail simulations. We argue that the ditterences 
are primarily attributable to resolution effects, specifically 
manifested in the distribution of supernova feedback en- 
ergy in different simulations. In high spatial resolution 
simulations, the added heat remains in dense regions and 
radiates away almost immediately, whereas in lower spatial 
resolution simulations a significant fraction of heat ends up 
in diffuse intergalactic gas where it cannot radiate away. 
Still, the differences in present-day WHIM fraction are not 
large between these two somewhat extreme cases for feed- 
back. Radiative cooling plays a minor role in the evolution 
of warm-hot gas because it is mostly at low densities, and is 
too diffuse to self-gravitate and self-shield into mini-halos. 
Other physical effects such as star formation, photoion- 
ization, and cosmic variance do not produce significant 
differences in the amount or properties of warm-hot gas in 
our simulations. 

Since the exact predicted fraction of WHIM gas is as yet 
sensitive to simulation details, it would be greatly benefi- 
cial to place observational constraints on this component. 
As mentioned before, this is a challenging task. S till, there 
are t antalizing h ints of detections of WHIM gas. Wang fc 
McCriy (1993) found an excess emission component in 
ROSAT data around T ~ 2 x 10^ K, which may b e arising 
from this diffuse component, ^oltan et al. (1996) detected 
an auto-correlation signature between the soft X-ray back- 
ground and galaxies, as would be expected if warm-hot gas 
was distributed in large-scale structures. Even more excit- 
ing is a possible direct detection of diffuse soft X-ray emis- 



sion associated with a filament of galaxies, by 3charf et al 



(2000|. Somewhat stronger evidence comes from a census 



diet the majority of WHIM gas is far away from galaxies 
and clusters, residing in the diffuse IGM. A promising av- 
enue to detect this more typical WHIM gas is via absorp- 
tion, as continuing observations with STIS aboard Hub- 
ble will detect many more O VI absorbers. Future X-ray 
satellites may be able to detect higher ionization absorbers 
such as O VII and O VIII that may also 1. race WHIM gas 
( [Hellsten, Gncdin fc Miralda-Escude 1998|) . The detection 
of this component is a key observational challenge, as the 
warm-hot intergalactic medium is rapidly become an inte- 
gral part of our understanding of the evolution of baryons 
in the universe. 
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